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ABSTRACT 


The growing number of large scale applications of network 
models and the availability of very fast solution codes make 
it attractive to formulate problems as networks whenever such 
models are adequate for the purpose. In this thesis, concep- 
tualization of, and notation used to express these models is 
based on the interpretation of physical flows of commodity 
through a network structure of nodes and arcs. As an aid to 
modelling, and to allow codes of varying specificity to be 
used, nine well-known Transformations are catalogued here for 
Gacy reference. 

Two recent results for special cases of the multicommodity 
flow problem are re-derived and in the case of (1) below, is 
Significantly extended: (1) The case with all capacitated 
mies in the network structure incident with one common node- 
(2) the case of a transportation structure with two sinks 
(or two sources). Using the network approach, these are 
shown to have equivalent network formulations. 

Lastly, a Transformation which uncapacitates a network is 
implemented in various ways into a contemporary solution code 


named GNET. 
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LIST OF SYMBOLS AND NOTATION 


eescet  OLenodes,pand directed ares, respectrvel ys 
which define the structure of a network. 


Flow aloweeare (1,1) from node 1 tounode j= 
Xs 1s taken as > 0, with flows entering and 
leaving a node having opposite signs. 


Upper bound on flow along are (1,}). 


Upper bound on sum of flows along a specified 
set of jointly capacitated arcs. 


ioe eselawer bounds, as Ua. sand Uso cues 
lippesenoundse. tre omitted, EL is understood to 
cme 1a): 


supply/demand at node 1. Used when it is not 
necessary to distinguish between a supply (T> 0) 
aidedemandstl < 0). 


AesuUppdy Of aunits available at node i- 

A demand of b units at node j, (b; > Oe 

DyesuUntt COSt Gk veransportation aleng sane sie 
Depending on comtext, signifies: 


a. Total number of Sources in the standard 
transportation problen, 


or 


b. lotal mumber of Sources which anresgeon- 
nected by permissible arcs to a particular 
node in a transshipment node. 


As for mM, a@ppilaed to sinks instead of sources: 
Signifies the commodity corresponding to that 


number in a multicommodity problem, e.g., x45 
mS tie tlowsoL commodity 2 along “are (154) 





ay 
L713 


Third subscript 


ale 


An example showing part of a trans- 
Shipment network. Note that the short 
arrows incident to the nedes indicate 
existence of other arcs and nodes not 
drawn--only the illustrated subnetwork 
1s-Of interest: Node jd asa source 
with supply aj and node 3 is a sink 
With, cdetand bz. Ares Mil. 2) sandals) 
have a jOint Upper capacity or eUeumints. 
and. (1,4) has as lower bound sOsn an. 
units. The curved lines intersect arcs 
to which the bounds apply. 


Signifies a member of a set, of multiple 
auwes Detween tWO MOdecemca 7, Xijk 1S 
the flow along the kth of a set 

Oimwares, 1rom node 1 togmodce wy. 





I. INTRODUCTION 


A. NETWORK MODELS 

A large and growing variety of problems in areas such as 
the transportation of goods, design of communications and 
mupeline systems, assignment of men to jobs, bid evaluation 
and production planning are adequately and accurately described 
by network models. 

These problems usually occur in the context of transporting 
a single homogeneous commodity from a set of points called 
'cource nodes" which generate supplies of the commodity sence 
factories, refineries) to another set of points called "demand 
nodes" with demands for the commodity (e.g., retails coresr 
end users), through a transportation system or "network" which 
consists of a set of notes (including the sources and sinks) 
connected by a set of directional routes or wanes (cvcess 
roads, pipelines). Each arc connecting a pair of nodes has 
a cost of transportation per unit of commodity, and perhaps 
an upper bound or "capacity" on the total flow of commodity 
along it. The objective is to find a minimum-cost flow pattern 
through the network which eatisties all demands. 

By network models is meant the special class of minimum 
cost (and maximum flow as a particular case) linear programming 
(L.P.) models of which the most general is the sSingle-commodity 


capacitated transshipment model, usually formulated as follows: 





Minimize ee OG... 3 X 
aS) 


Subject to ? he ahs qT, HOl Caen motenia Ehis 


equation expressing the 
algebraic relation between 
incoming and outgoing flow 


at the node, 


and [Geer << Ue £0Y Cach anew 
1) 1) 1) 


Non-negativity of the sae and equality of total supply to 
total demand is necessary in the above equality formulation 
Siemtne model. Positive lower bounds, ee on 2s are easily 
dealt with by simple arithmetical transformations to yield 
Mertables with zero-valued lower bounds. 

Other models in this class are more specific cases of the 
above formulation, and include the capacitated and uncapacitated 
transportation model and the personnel assignment model. All 
these models are well-known and widely described, for example, 
meker. [| 9]. 

Although the models are conceptually and algebraically 
samole, the real-life problems they describe are typically 
huge, and it has been found necessary to find economical means 
of solving problems with, for example, 50,000 nodes and 500,000 


arcs. 


B. VIEWPOINTS FOR NETWORK MODELS 
The notation used in expressing network models is important 


because it influences the development and implementation of 


10 





network algorithms, and is relevant to the ease with which 
Gereain theoretical results are conceived and proven. 

One way of viewing a network model is as a straightforward 
meee, With slacks introduced to convert the inequality capacity 
@emstraints to equalities, the constraint expressions in the 
above formulation may be economically written in the usual 
Meerix form MX=-T. 4L.P. solution codes accept input and manip- 
ulate data in matrix and vector form, and theoretical results 
arrived at using this approach involve operations borrowed 
from linear and matrix algebra. 

ANOTHeET Viewpoint, which may be called the network approach, 
1s derived from the actual physical nature of many problems-- 
that of finding an optimal pattern of flows His through a 
network structure consisting of directed, capacitated (perhaps) 
arcs (i,j) which connect nodes at which are demands/supplies 
T;- hits Viewpoint deals witheares Gather than vectors. and 
interprets the L.P. constraint equations as expressions of 
flow conservation at a node. 

The network approach is more natural than the use of the 
notation in the above formulation of a network model, since: 

PP onnpcacti cal Network problems, 1t 1S. tanely tHe 
case that every node pair is joined by an arc. There is also 
mires possibility of multiple arcs between a node pair. It is 
thus more appropriate to formulate the problem in terms of 
the set A of arcs which are actually present in the network 


Seructlture: 


1] 





° ° ° a 
Minimize (i,j)eAcij* a3 
; Sarr ease eo ee a jie 
subject to Ree en 1 fal (Ce eke T., for each 


node i, 


and Perens UU... L£Or seach (ieee ek. 
13 13 13 


b. Data input and manipulation in all efficient 
network codes is in arc form. The storage and manipulation 
Gmeeraree Sparse Matrices implicit in L.P. codes is by compar- 
ison inefficient and cumbersome. 

In view of the growing size of real-life applications, a 
principal practical result of the network approach has been 
the development of solution codes which are much faster than 
ime codes for problems of equal size. Although they follow 
pmeeeisely the SIMPLEX algorithm as do L.P. codes, network- 
besed Codes accept, store and manipulate data in a manner 
which is based on the interpretation of a basic set of vari- 
ommes as a tree which spans the nodes of the problem (Koopmans 
[6})}. "Tree" is a term from graph theory which describes a 
network which has only one (unique) path between any two nodes, 
disregarding arc orientation. 

The differences between the two viewpoints lie in visual 
conceptualization, terminology, and solution implementation, 
rather than in any difference in mathematical foundation. 

The solution algorithm is the same in both, and as an example 
of their commonality in mathematical basis (which 1ies in 


linear algebra), it is noted that the network transformations 
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introduced below may be viewed simply as linear operations 

on the matrix M and vector T. In many respects, the two view- 
points are complementary in rigor and ease of intuitive appre- 
SEdcion. Ihe physical interpretation of the problem in the 
network approach does make it particularly acceptable to 


non-analysts. 


C. THE MULTICOMMODITY FLOW PROBLEM 

The general multicommodity flow problem is a minimum-cost 
transshipment problem in which a number of distinguishable 
commodities flow along the capacitated directed arcs of a 
network structure in accordance with separate demands and 
supplies for each commodity at the nodes, and an arc's capacity 
apetaes to the sum of flows of the commodities along it. There 
Meowpcen recent interest in finding special cases which may be 
transformed into network models in order to take advantage of 
fast network codes, the more so because a problem with r conm- 
modities has approximately r times the number of variables of 
aeoingle commodity problem with the same network structure. 
D. IMPORTANCE OF NETWORK TRANSFORMATIONS 

AND CONTRIBUTION OF THIS THESIS 

Because of the size of many practical applications, and 
the availability of contemporary network-based codes which 
@amepe faster than L.P. codes by a factor of 100 or more, it 
is obviously advantageous to formulate a problem as a network 
model wherever possible, 

ittawitt be noted tat. GCapacitated network models have 


Capacities which are imposed only on individual arcs. 
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There are many problems which have capacities imposed on nodes, 
ana@eon sets Of arcs (a “joint capacity" being an upper bound 
on the sum of flows through a specified set of arcs). Such 
problems cannot directly take advantage of network-based codes, 
omaemust be re-formulated or transformed into an equivalent 
Capacitated or uncapacitated network model (there are, however, 
@-feain types of joint capacity described later which do not 
Meamit this transformation). 

Ht as also important to be able to transform network models 
amemast themselves (é.¢g., transshipment to transportation, 
Capacitated to uncapacitated) to suit the specificity of a 
particular solution code, since not all network-based codes 
eam sSOlve the general capacitated transshipment problem. 

For the above purposes, there have been a series of trans- 
moamations proposed and used over the past 20 years or so to 
deal with networks which are jointly capacitated in a number 
Of ways. The main work of this thesis lies in: 

1. Cataloging these transformations into a single refer- 
ence as an aid to modelling. 

2. Developing a unified approach to the many transfor- 
feaereons by Showing that they may be constructed by selective 
sequential application of a very few basic transformations. 

Be oRe-derilVvation Of two recent results (and in one ease, 
Saoenificant extension) involving special cases of the multi- 
commodity flow problem, using the network transformation 


maeproacn. Ihe two cases of interest here are: 
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a. The two-commodity transportation problem where 
all capacitated arcs are incident with a single node (Rebman 
Meer. Here, the original result is significantly extended 
to the r-commodity transshipment case, in which the permissible 
mimes need not be identical for all commodities, and an arc's 
@am@acity need not apply to all commodities flowing along it. 

b. The multicommodity transportation problem with 
emener two sources or two sinks. 

ties sapplaecat lon 15) an 1RUStratron- on the wlilveysot 
meaework transformations in constructive proofs of the existence 
Gea equivalent network model for a given flow problem. 

4. The implementation in various ways of a transformation 
designed to Uncapacitate a transshipment network, using an 
existing network-based code GNET, which is capable of solving 
@mewcapacitated transshipment problem. The speed of solution 
with the transformation implemented is examined to see if the 
Sawanie in required memory space for a projected uncapacitated 
version of GNET would be worthwhile as an offset against any 
reduction in speed which may be caused by implementation of 
Maem Ganstormation. 

It is emphasized that the transformations derived here are 
not new, and most of them have their origins buried deep in 
Mme tOlkore’’ of this field. This thesis will attempt to 
give credit where it is known, but an apology in advance to 


Originators not mentioned herein is not out of order. 


aie. 





E. SOME IMPORTANT PROPERTIES OF NETWORK MODELS 

iiewconstraint matrix M in the formulation MX=T of a 
network model has the property of total unimodularity, for 
Waren a necessary and sufficient condition is that every 
Square nonsingular submatrix of M has a determinant of value 
mone. tIhere other kinds of models which have this property, 
but all practical totally unimodular problems have so far been 
shown to have network model formulations. 

fron the L.P. viewpoint, the transformations derived in 
this thesis are viewed as reformulations of the original 
problems in such a way that the matrix M of the transformed 
problem becomes (or remains) totally unimodular. 

memajor consequence of total unimodularity is that eweny 
Deeete SOlUtIOnN (and an optimal solution) is integer valued, 
8ivyen integer T in MX=T. This is of importance in many real 
m@easioeems fOr which only integer answers make sense. It is 
pueeee tive that if the optimal solution to an L.P. problem is 
unique and fractional, then the problem cannot be formulated 
as a network model. 

Miere are other classes of L.P. problems which are inter- 
pretable as finding optimal flows through a network structure 
but which are not totally unimodular and thus cannot be formu- 
ited aS a network model as defined here, except in special 
cases. Examples are the network flow problem with gains, and 
the general multicommodity problem, which can in fact be thought 
of as a single-commodity problem with an expanded network struc- 


ture made up of disjoint sub-networks (one for each commodity) 
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identical to the original structure, with a certain type of 
joint capacity which does not permit a network model formula- 
tion. Such flow problems are not termed "network models" in 
the sense used herein. 

A network structure is said to be bipartite if its nodes 
can be exhaustively divided into two disjoint subsets such 
that every arc has one incident node in one subset and the 
Sener incident node in the other subset, and all arcs (directed) 
are oriented in the same direction from one of the subsets to 
@eesother. The structures of the transportation and assignment 
Peoutwems are bipartite. 

In the terminology used here, nodes through which the com- 
femmnty 1S permitted to flow en route to other notes are called 
transshipment nodes. A node at which total incoming commodity 
is constrained to equal total outgoing commodity (i.e., at 
which demand = supply = 0) is called a "pure transshipment" 
Meee. oOuUrces Or Sinks which are not also transshipment nodes 
eemediied “pure source’ or "pure sink" nodes, otherwise "demand 
Paaiosnipment'’ or “supply transshipment" is used. When it is 
unnecessary to make distinctions, the word "pure" will be 


omitted. 
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Ii. NETWORK TRANSFORMATIONS 


Of the following Transformations, the first three may be 
regarded as basic in the sense that they are used in the con- 
Struction of almost all the succeeding (more complex) trans- 
emma tions. 

Transformations are named and referred to by Roman numerals, 
pmemeare spelled with a capital “T" to distinguish them from 


normal usage of the word, in this and subsequent sections. 


A. USES OF TRANSFORMATIONS 
These Transformations have been constructed with one or 
more of the following uses in mind: 
“aero Simplify a complex capacitated structure (e.g., 
transforming a group of arcs with a joint capacity 
mito arcs with individual capacities). 


meee Oo uUnCcapacitate a capacitated network. 


Gee l.O transform a transshipment network to a transportation 
me twork. 


d. A combination of some or all of a, b, and c, above. 


e. To simplify a network, in the sense of reducing the 
number of arcs and/or nodes. 


The broader purposes of using transformations lie in: 

1. Transforming a problem into a form suitable for use 
Tieethe particular type of network-based code available. 
The specificity of application of a code may vary from the 
Mmeapacitated transportation problem to the general trans- 


Shipment problem. 


18 





For problems of equal size (i1.e., with the same number 

Of arcs and nodes), a more specific code should be faster than 
aeeeticra) code. On the other hand, transforming a network 
into a form which can be handled by a more specific code 
almost invariably carries the penalty of increasing the size 
of the network due to the addition of artificial arcs and/or 
nodes. Transformations which reduce the size of a network 
Carry the reverse penalty of increasing the generality of the 
moe nrequired, C€lther because they introduce more complex 
Capacitating mechanisms or because they transform a transpor- 
maewen Structure to a transshipment Structure. There is there- 
fore an inverse relationship between the size and the generality 
of the many forms a problem may take. 

Z. Use in the constructive derivation of results, as 15 


@eme later in this thesis. 


Be EQUIVALENCE OF NETWORKS AND INVERSE TRANSFORMATIONS 

The idea of "equivalence" in transforming one network into 
paomener JS detined in terms of being able to obtain the optimal 
solution to the original network by performing post-optimal 
iemenmetic Operations on the optimal solution (objective func- 
tion and/or flows) to the transformed network. It is necessary 
for equivalence of two networks that their objective functions 
darecer Only by a known additive or multiplicative factor. 

Equivalence also implies "invertibility" of a transforma- 
mion tO Obtain thé original network. However, given only the 
transformed Metwork, the applicability of the inverse trans- 


neEmimabtOn 1S cE Lieult to See an most cases. .In face; it is 
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rare to find a network with the precise features necessary 


for application of many inverses (e.g., zero-cost arcs). 


SeeeeeEXtolENCE OF EQUIVALENT NETWORK MODELS 

In order for a network flow problem with joint capacities 
to have an equivalent network formulation, it is necessary 
for its joint capacities to satisfy the following well-known 
@omaations (e.g., Dantzig [| 3 J): 


Condition I. At any node, any two joint capacities imposed 
Om ats incident arcs either have no arcs in common, or the 
arcs of one joint capacity are a subset of the other joint 
@upacity (1.€., any two joint capacities at a node apply 
Mmommare SeCLS which are either disjoint or such that one are 
See 1s a subset ("nested") of the other). 


awrclon f). A joint Capacity is not imposed on arcs any 
mgomor more Or which are completely node disjoint (i.ée., 
mount Capacity should be applied to a set of arcs which 
has at least one common incident node). 


Menaition Ili. A joint capacity at a node is not imposed 
on a set of arcs which has some member(s) oriented in the 
opposite direction to some other member(s), with respect 
meme node. Ihe exception is when all the arcs are inci- 
dent between the same node pair. In this case, only one 
Ot the arcs will have a non-zero flow in any basis. Thus, 
the problem can be reformulated with each arc individually 
Capacitated with the original joint capacity. 


more | illustrates violation of Conditions I and il, 


megmeeroure 2 illustrates violation of Condition II1. 


Uy 


Ptge. flGwc 
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When a joint capacity appears visually to violate Condition 
II, it is sometimes profitable to look at both ends of the arcs 
iiedmestion. For example, a first glance at Fig. 3 indicates 
emesis) ViOlation of Condition II, but when the capacity is 
Oman aS in Fig. 4, it is seen that this is not the case. If 
Pree o had a more complex structure, its equivalence to Fig. 4 


iyeemot be at all clear. 


e O O 


{oN \ 7 
Pe / s oe / I) 


MiG ¢ Fig. 4 


This approach can be particularly helpful in problems with 
Supmet works which have transportation structures, where the 
known sum of flows at any node can considerably simplify a 
women capacity. For example, consider the equivalent capa- 
Cities in the identical networks in Figs. 5 and 6. Note that 
Fig. S appears to violate Condition II. Fig. 6 is obtained 
bpaene ting that in Fig. 5S, U includes all arcs incident with 
the pure demand nodes 3 and 4. Since the flows along these 
Mees must tOtal exactly bz cs Dy» then the remaining arc 
(1,5) covered by U must then have a capacity U- (b, + b,): 
21 








Fig. 9 | Fig. 6 


In the rest of the thesis, a joint capacity in any problem 
is assumed to satisfy conditions which ensure a feasible solu- 
tion, and also conditions which ensure that the Capacity 1s 


Meme redundant (1.e., binding). 


D. THE TRANSFORMATIONS 
ieee [ranstormation I 
This Transformation reverses the direction of a capa- 
ereated arc. 
Suppose (1,2) is an arc ina transshipment problem, 
Teeetinit cost of flow c, and such that Xi Us Eaen toa 


nodes 1 and 2 may have demand or supply T,, I>. See Fig. 7. 





Aids. 7 


Ze 





Cometderwcie eCauivalent network in Fig. 8, which can 
be regarded as the result of introducing a counter-flow of 
magnitude U in the arc, which gives a resultant flow of X51 


eee X Mine maIrectTon opposite to the original flow x 


eZ 


and making the cost negative. 


DEP 





Pidesc 


Algebraically, this is simply equivalent to substitut- 


miei - xX ii place Of x Wherever it oCeumsS inthe alow 


‘eal 12 


conservation equations and objective function of the original 


problem. 


Original Problem Transformed Problem 
Node 1 flow oO) Oe 
servation: se “Xy5 + T) = 0 =e -(U - X54) + T) = 0 
Or ... + X54 ba - U) = 0 
Node 2 flow con- 
servation: see tXyo + T, = 0 Pres ne X54) a T, = 0 
OF «4. “X5oq 2G + U) = 0 
Pmemcapacity: OCa X15 < U OS Xo, U 
On US oe 0 
Sieyective function:... CX,5 t oe Pe X5,)+ 


7. 





Tie soute@en toethe Orleinal problem is obtained from 
the solution to the transformed problem as follows: 

a. X15 * = Xo4° Other flows are identical between 
the two problems. 

DeeeradmeclictO “themoujyective Of the transformed prob bene 

lGansmomiarvone tis .0r little interest by i1tselt. bwe 
memasectul as part of a more complex transformation. 

2. Transformation Ila and IIb 

GiicimamomiercCOuSubaINt, OVEr a SUBSEL Of Ehevares 
incident to a node, this Transformation (actually a sub-class 
wee nanstormations) transfers the capacity to a single (arti- 
feed) zero-cost arc, leaving all other arcs uncapacitated. 

wena oldu, tOlating Condition 11), the jointweons trai 
mist be imposed on the algebraic sum of flows in the arcs. 
Pmews MUSt therefore have attached signs appropriate to their 
fect atlOn with respect to the incident node. 

Constraints having sums which include flows with oppo- 
Site signs mist be clearly specified as to whether they are 
permitted to have negative lower bounds. 

Cons1cer tne transsminment Network in Fig. 95 wither 
mommesconstraint as shown. Note that a constraint which in- 
4 


eiuidees the sum xX + for example, is not permitted 


Ze” See 
momeene reason given above. The equivalent network is shown 


terse. 1. 
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with constraint: 
0 £ X94 ee en) 
Fig. 9 Fig. 10 
A positive lower bound is also permitted, with or without an 
upper bound. | 
id X5, + Xz, 7 Xqy iS NOt Constrained to be non-negapave, 
heel results. The artificial arc pair is necessary to allow 


the sum X54 “ Xz, 7 Xa 4 freedom to flow in either direction 


(feeees according to its sign) between nodes 1] and A. 





Ze 





PecOnsctraimt et the £orm typaried by X57 * Xz, 7 Xqy 
emve when written as Xo) + X3y S Xyq> Can Degeas ily  scen to 
be useful, for example in a road network where the total in- 
coming traffic along some subset of roads at a busy town is 
moaned to be no greater than total outgoing traffic along 
some other subset of roads, while the more general form 


“3 - Xyq > U would require a positive difference between 


oor” 31 
the two sets of traffic flows. 

Transformation IIb is a special case of IIa, and is 
commonly known as the "capacitated node." It arises when an 
upper (and/or lower) bound is imposed on the amount of flow 
@aesane through a node, which, in the case of a pure trans- 
shipment node, is equivalent to putting a joint capacity on 
mueeene arcs entering (or leaving) the node. Transformation 
memeean then be applied directly. For a source or sink node, 
the demand or supply at the node may be included as part of 
the total flow passing through. Figs. 12 and 13 show the 
Transformation for a source node with an upper bound U, which 
maiemaes the supply as part of the total. If the supply is 


Hemeeo pe included, then the supply ay 1S placed? atenoders 


instead of node A. 





Pied Fig. 13 
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A set of jointly capacitated arcs at a Modes 1s a 
Special case of a jointly constrained set, when all arcs in 
the constrained set are oriented in the same direction. 

3. Transformations IIla and IIIb 

Given a capacitated arc, these Transformations consist 
of Transformations II and I applied sequentially to create an 
uncapacitated transshipment structure with one additional 
eeiciicial aces aid med ec. 

Consider the capacitated arc (1,2) ina transshipment 


network, Fig. 14. Applying Transformation Li gives VE re weer 





rigs 12 Fig. 15 


Applying Transformation I to reverse the. are. Clas 
Paveo Fig. 16. It is seen immediately from flow conservation 
aienede A that flows Xay = U - Xq5 and Xn = X42 GCannot execea 
U. Thus, the redundant capacity on arc (A,1) can be removed, 


with the end result in Fig. 17. Call this Transformatiom, [iia 





Fig. 16 Fig. 17 (Transformation IIa) 


aa 





No modification to the objective function is necessary 
femobtain the solution to the original network. The original 


fliow x VomEOUNGdmaEneectly an arc (A,2Z). 


12 
Alternatively, the original network may be transformed 
memenat shown in Fig. 18, instead of Fig. 15. Reversal of arc 


(A,2) with Transformation I, and using the same reasoning as 


apove, gives the end result in Fig. 19, which has reversed 


@aesecompared to Fig. 1/7, and the original flow X12 is found 


Gapeeetly in arc (1,A). 





Fig. 18 Fig. 19 (Transformation IIIb) 


Peelcomn ceed eines the twomend Tesults typatied= Imeem 
17 and 19 may be transportation structures, if the demand/supply 
magne Orientations at nodes 1 and 2 are appropriate, since 
node A is already either a pure source or sink node in the 
transformed networks. 

wemelransformation IV 

This Transformation, applied to a jointly capacitated 
elmsct Of arcs at a node in a Capacitated transportation net- 
work, produces an uncapacitated transportation network with 


mmemmaddition of two artificial nodes and arcs. 
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Suppose we have part of a transportation network as 
in Fig. 20. Application of Transformation II yields the 
equivalent transshipment network in Fig. 21. 


“by 





Fig, 2D» Fig. dl 


itenstormation 1flb applied to arc (1,A) gives the 
network in Fig. 22, which can be seen to be part of a trans- 


WeieeatiOn network when re-arranged as in Fig. 23. 





Fig. 22 Fig. 23 
Nowmodif£iCatien LOwthe Objective function isonece. sam, 
and the original flows Xy9° Xz are found’ an ares (AyZ eee eo 
mempectively. 
Suppose the original problem had a jointly capacitated 


SeenOt arcs entering a sink. Transformation Il, and then 
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Transformation IIIa (instead of IIIb), would produce a 
Simeteror image’ to Fig. 23 above. 

This Transformation is of course applicable to indi- 
imieually capacitated arcs. If ail the arcs incident at a node 
ieee a transshipment network have a joint capacity U, then for 
iewomoe Dinding, it 1S required that U < a.. This case makes 
Bemse Only if it occurs in a network with total supply Za, > 
total demand oe Otherwise U < a. makes the problem infeasible. 
Additionally, only a source node may be capacitated in this way, 
since a binding One canine U < Ds at a sink immediately causes 
infeasibility (the case where Ub. > bas is not considered). 
Mmetemis Special case (Fig. 24), a, may be replaced cirecet 
Mameeco £1ive the uncapacitated Fig. 25, which requires no 
fae icial arcs or nodes. However, in the standard procedure 
for the case of ha. > Pleas where all the sources are connected 
at zero cost to a dummy sink with demand -(Za, ji He this 
itwemstormation for the special case in Fig. 24 necessitates 


Suse cution of U for a 


} in the demand - (4a, - PDs): 


U 
ay U 
Fig. 24 Fig. 25 
me liranstormation V 


A network with a node at which one joint capacity 1s 


Pmesteda’ within another in the manner described in Condition I 
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in Section II.C, is transformed to a transshipment network 
which includes two individually capacitated arcs. In this form 
it may be used directly with suitable codes. If the original 
network was a transportation structure, the transformation 
may be carried further to give an uncapacitated transportation 
network. 

Suppose we have part of a transshipment network as in 
Fig. 26. Necessarily, U, 7 U, Od: feasibility and non- 
redundancy of U,- Application of Transformation II to arcs 
capacitated by U, > and then again to arcs capacitated by U5; 


pews Fig. 2/7. 





mide 26 Fig. 27 


Suppose now that the original problem was part of a 
Mmemiscportation network, and it is desired to obtain an un- 
Capacitated transportation network. Beginning with Fig. 28 
maimeneis analogous to Fig. 27, application of Transformation 
ieee o arcs (1,A) and (A,B) eives Fig. 29, which can be 


mewamecad in the usual bipartite form as in Fig. 30. 


oo 








Fig. 28 


Fig. 30 


ee 





Mimeitilctmiile. = 'c/ Or rip. S50, no Change to the obj ec- 
tive function is necessary, and any original flows ae are 
found in that arc which enters node j in the transformed 
Ge blem. 

6. Transformation VI 

This Transformation is one means of formulating a 
transshipment problem as a transportation problem. It is 
found, for example, in Wagner [9]. 

In this transformation, every node in the transshipment 
problem which has all its incident arcs either entering or 
Memvine (1.€., 1S either a pure sink or source node) is left 
unchanged. The following transformation is then applied to 
every transshipment node, a typical example of which is in 
memes, ithe node A can be split into two nodes, one with ail 
the original arcs leaving (node Al), and the other with all 
Ommminal arcs entering (node AZ). A zero-cost arc (AZ, Al) 


mmteaoed aS in Fig. 32, which is clearly equivalent to Fig. 31. 





rg... Ol Fig. 32 


in anticipation of using Transformation I, a super- 
mretis Capacity is imposed on arc (AZ, Al). A safe Capacity 
Pemuse 1S Obviously the sum of all supplies in the problem 


for any larger number); since no sum of flows in the set of 
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arcs entering (or leaving) any node can exceed this number, 
which is denoted by B, and given the interpretation of a 
Sbmerer stock at node A" in Wagner [9 ]. 

Applying Transformation I to reverse arc (A2, Al), 
Peeeemence Obtain a bipartite structure, we obtain Fig. 33. 
Paemoae AZ, the flow conservation equation is now x NIA? = 
B - iZA1* i, A2? which ensures XAq,A2 <b) The conse cain cao 
arc (Al, A2) is therefore superfluous, and is removed to give 
finally Fig. 34, in which Al and A2 are pure source and sink 
Modesmrespectively. Application of this transformation to 
ey transshipment node, together with the unchanged pure 
Bemeecrand Sink nodes in the original problem, will clearly 
Petanan cquivalént transportation problem in which the original 
miowemare directly identifiable and the objective function 


unaffected by the introduction of zero-cost artificial arcs 


such as (Al, A2). 





FG .ao3 Fig. 34 


mire above transformation is not affected by capacities, 
fourveeor individual, on the arcs in the original problem. 
Capacities in the original problem are imposed on the corres- 


momaing arcs in the transformed problem. 
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7. Transformation VII 

This Transformation is applicable to a pure source 
node which delivers flow to precisely two other nodes. It 
Ewamimates the source node and its incident arcs, at the cost 
Sememtroducing an artificial capacitated arc. It can there- 
fore be used to reduce the number of nodes and arcs in problems 
fieneethe special structure mentioned above, if a code capable 
of handling capacitated transshipment problems is available. 

BUPppOse the Network ImeFileg. $5 1S part Of a trans- 
Sampment problem. This is re-drawn in Fig. 36, with a super- 


miMous capacity a, imposed on arc (1,2). 


1 





Fig. 36 
mppiyine Transformation £ to reverse arc (1,2) yuellds 


fiers which is algebraically simplified in Fig. 38. 







a 2g: mG 
-b 
ta -b Z -b 
i] 3 +q ~C € 3 
ay 1 12 4, ne 
mlela cy Fig. 38 


Node 1 in Fig. 38 is clearly superfluous, and its 
removal yields finally Fig. 39 which is equivalent to the 


Gurcinal problem provided: 
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a. X)2 in the original problem is obtained from X52 


in the transformed problem. 

b. X19 in the original problem is obtained by sub- 
tracting X52 from ay: 

fae inesoricinal Opyjective function contained casts 


C1 2%13 if C1 9%4° The transformed objective function contains 
pee 12) ~~ 9 13%a3) “a 2%13 ~ ©135%13 ~ ©1242 ~ *q2)! = 


7 C thas the original cost 1s obtaaned 


iets = ~12°12 ~ “1~12° 
by adding a)C1> to the objective function of the transformed 


problem. 





Fig. 39 


This transformation can be interpreted as if all of 
supply a, were delavered firstuseeuwnode 2, %aty cost a,Cy>° Sub- 
Peomemtly, each unit increase in X52 (= X13 in the original 
problem) will increase cost by Cy2 units and decrease cost by 
C15 Upees. This Can be seen in Fig. 3 to be precisely the 
Samemmechanism for a unit increase in X13 in the original 
problem, and arises because flow conservation at node 1] forces 
feeminc decrease in X19 LOm@wcVery Ulit TiChease 1n X13" 

weet ranstormation VIII 

has Transhermationm is applied to effectively remove 

[7-e€., reduce to zero) a lower bound on the flow along an arc. 


Suppose we have arc (1,2) in Fig. 40 with both an 


Mpper and lower bound. Since at least L units must be delivered 
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along the arc, we can simply apply a tranformation which 
assumes that L units have already been delivered, leaving 

the balance X19 7 L (> 0) to be determined by the solution 
algorithm. The upper bound has also to be correspondingly 
modified as shown in Fig. 41, where X19 relers to flow molec 


transformed network. 





Fig. 40 - Fig. 41 


ihe original flow X15 1s obtained from x9 + Leama 


the original objective function value by adding Lc Comene 


eZ 
optimal solution of the transformed network. 

In the extreme case when U = L, equivalent to specify- 
ing Xo = U, we can eliminate arc (2) tO mob tain noe 


euicabie operations must of course be carried out to obtain 


PiemOoraeinal solution. 


Fig. 42 When U=L 


ome transformation 1X 
This transformation is applicable to a pure source node 
which delivers flow to two nodes via capacitated arcs. It 
reduces the number of Capacitated arcs to one, and thus may be 


peerul in reducing the complexity of a given problem. 
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Suppose Fig. 43 is part of a transshipment network. 
We have X73 cS Uj, 3<=> 4, - X)3 2 a, - U1 3<=>X1 5 > ay . Uj: 
Maewarc (1,3) can therefore be uncapacitated by imposing a 


emer bound on flow in arc (1,2), resulting in Fig. 44. 








as 
24427 
a3 a 


ay 





Fig. 43 | Fig. 44 


Eeplacation of Transformation IX results in Fig. 45 with 
Sumeerone Capacitated arc. The original flow X19 is obtained 
OM X19 + a, -U,3, and the original objective function value 


by adding (a, - Uj, 3)¢)> post-optimally in the transformed problem. 





Fig. 45 
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Pel eewNO SPECIAL (CASES ©OF THE 
MULTICOMMODITY FLOW PROBLEM 


A. GENERAL DISCUSSION 

When viewed as an extension to the capacitated single- 
commodity transshipment model, the multicommodity model has 
naturally wider applicability in more complex and realistic 
Sieeuations. Apart from the more obvious applications of 
Shipping different goods through a distribution system, it 
1S sometimes necessary to keep track of groups of items which 
teweeoe physically identical. For example, some traffic 
asomemment problems require that the identity of travellers 
going between various node pairs must be kept. The model is 
Mocamin problems which require finding optimal flow patterns 
(quantity of flows and routes) for each of a number of dis- 
memeuishable commodities in a capacitated network structure. 
iWeweriteria of optimality is minimum total cost, which in- 
cludes the maximum sum of flows criteria as a special case. 
Pecapacity On an arc takes the form of an upper bound on the 
sum of flows of all commodities along that arc. Only directed 
ieseare considered in the following results, since a capaci- 
eeeeae undirected arc in a multicommodity problem does not in 
general have an equivalent network formulation because, even if 
the usual single-commodity formulation of two jointly-capacitated 
fumes Oriented in Opposite directions between the two incident 
nodes is resorted to, there will generally be flow (of differ- 


Cimecoumnodi1ti1eS) in both ares, unlike in the single commodity case: 


ou 





For any one of the commodities, the flow in only one of the 
iPvemares Can be at a non-zero level. 
The general minimum cost multicommodity transshipment 


peeprem with r commodities is usually formulated as follows: 


. ae 2 Roe ok 
Minimize VE ora) ED, © area 
kel ij #1) 1) 
Subject to: 
Tx. = rx - TX 5 Je) eclel Euler 
j.7J fi oak i) 


toe kk ne 
pees sae = Vis Menor Cache. j., 


and es > oO fOr geaCli days jigs: Ke 
13 
Additionally, it is assumed that ies = 0 for each k, 
Hl 


The general multicommodity problem above is non-unimodular 
(thus, even with integer e Ui? tne solution will inieemenall 
be fractional) and is solved as a linear programming problem. 
For large problems, therefore, the comparative slowness of 
L.P. techniques is again a limiting factor, more so because 
an r-commodity problem is approximately r times the size of 
a Single-commodity problem with the same node-arc structure. 

In each of the following two subsections is provided (as 
Part of a constructive derivation) a transformation of a 
special case of the multicommodity problem into an equivalent 
Single-commodity capacitated transshipment problem. Each of 
these two transformed networks (or with further transformation 
into equivalent uncapacitated transshipment/transportation 


networks) may then be solved directly and advantageously with 


network-based solution codes, yielding integer solutions. 
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epecitically, the two results are?:- 
a. A network-based derivation and significant extension 
of a result by Rebman [7 ] which originally considered the 
2-commodity transportation problem where all capacitated 
arcs are incident with a single node in the network struc- 
more , 
b. A network-based simplified derivation of a result by 
Evans, Jarvis and Duke [4], which considered the r-com- 
modity transportation problem with 2 sinks (or 2 sources). 
ime Original results were obtained by using matrix theon., 
It is believed that the derivations below will show elleaqsiy, 
the advantages of the network approach with respect to sim- 
paememty, Visual clarity, and greater appeal to intuition. 


They also illustrate the use of network transformations as 


part of what amounts to a proof of unimodularity. 


fee oULT I-=-THE “BUSY NODE*” PROBLEM 
Rebman [7] showed that an equivalent single-commodity 
Memwwoerk formulation exists for the 2-commodity transportation 
peebpbem in which the capacitated arcs from a tree with at most 
Simeminterior node (i.e., all capacitated arcs in the problem 
are incident with a common node--the “busy node," which is 
eeener a source or sink node). 
Tt will be shown here that his result can be extended 1n 
the following ways:- 
weer 2 2 commodities. 
b. The network may be a transshipment network, in which 
tine “busy node" may be any node (source, sink, Or pure 
transshipment). 


See ttiltiple arcs between nodes are permitted. 


d. Permissible arcs need not be identical for each com- 
modity. A permissible arc for the kth commodity is defined 


4] 





as an arc along which the commodity is permitted to flow 
mimomposttive Level at a finite cost. 


e. Joint capacities and lower bounds may be imposed on 
mise ts Of the set of arcs incident with the busy node, 
Maovwaded the three conditions in Section I1.C. for an 
equivalent network to exist are satisfied. This includes 
the requirement that each member of a capacitated subset 
of arcs be oriented in the same direction with respect to 
the busy node, for every commodity. However, disjoint 
@epdeiitated subsets of arcs may be oriented in opposite 
directions. 


meet 1S not necessary that a capacity on an arc applies 
Bomene sum of flows of all commodities along it--1.e., a 
Sryoreity on the sum of flows of some subset of the r com- 
MomeericsS 1S permitted, the rest of the commodities either 
being delivered without bound along the arc, or having other 
joint capacities imposed on nested or disjoint subsets of 
commodities. 

Mirere 1s of course no requirement that all arcs incident 
fameeeac bUSy node be capacitated. It 1s required that equal- 
ity of total supply to total demand holds for each commodity. 

The result has obvious application to multicommodity 
erarsohipment problems with one particular central terminal 
Wikemerequires incident routes to be capacitated in order to 


Metdmaze congestion. Pictorially, a typical busy node coulda 


appear as in Fig. 46. 


Bee! ' Sil r 
b35a3 | Apres rdy 
X 1 is 
BY 7% ! 
Bb! bh 2 2 aime 15) 





are 46 


42 





The busy node B shown above is a source transshipment 
node. the heavy-lined arcs are capacitated, with lines drawn 
across groups of arcs which have joint capacities. 

Algebraically, the problem is formulated below, using 


the following notation: 


N The set of nodes of the network. 

A The set of arcs of the network, defined as 
node pairs. 

aX Gye thtesset Of permissible arcs for the kth 
commodity. 

A(i) Greene Set Occarcs incident to node 1 

A(1,j) GA, the set of multiple arcs from node i to 
mode) i. 

A[B] STi etne secu. Of Capacitated adres ineidenmtmad 
busy node Bb: 

uae, 13%) © A(i,j), the gth multiple arc from node i to 
OES ye 

kk | ' | 
Os 40 ane Unit cost and flow of kt" commodity along arc 


(Gisiae)) 2 (AK lite) )). If the third subeeruer 
is absent, it is understood that there is only 
Oieware (i. ) 2 


ies problem 1S: 


k _k 


— 2 
ijn SUCm that CF ie *aga1? 


iG 
Minimize es 
(le eee elie 9.) 


Ke L 


Su yeet FO: 
k — 
Xije OP debris jo “lee, Ky 
: x 2 Kk 
211,532) eAK Mai) 19% ninf(njismeakaaciy ™™ © 


all 
= Oe ee ele Kae 
i 


and at node B, 
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> ce Be 


x* 
gist ye 
kI (1,532) © AK ABI 


each (i,j;2) € A[B] 


feeavoid further notational complexity, the possibility 
Smepoimt Capacities or capacities which apply to only some 
commodities along an arc are not reflected in the above 
memm@eation. These extensions are dealt with later. 

dlc Constructive Detivattou Or the Result 

Suppose r = 3 for a transshipment network, part of 

hen 1S Shown in Fig. 47. Busy node B is in this example 
a sink transshipment node with demands bp. 


those nodes directly joined to B are shown, and in fact only 


POR Clarity jae, 


these nodes and the arcs joining them to B are visually rele- 
womemeorthne proot. Individually capacitated arcs are used 
Memo lacity. We consider here the case where permissible 
arcs are identical for each commodity, and the capacities apply 
to the sum of flows of all commodities along an arc. [It is 
seen later that the proof is unaffected by these simplifying 


assumptions. 


B3 
Vip 
U 1 2 3 
2B “bp »-bp »-bp 
ema? 
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The network is equivalent to r = 3 separate networks, 
one for each commodity, and each identical in structure to 
the original network, with capacities applied to appropriate 
Stieceacross the r networks. Let node aS in the network 
Comresponding to the kth CONMOGITTY , —COrrespold £0, Node M104 
the original problem. Fig. 48 results. We now have essentially 


identical subnetworks coupled by capacity constraints. 








commodity 1 
with capacities: 
3 
k 
pee oe eee 
k=] 1B 1B 
commodity 2 2 ae VoR > 
k 
YX ene 
aa B3 B3 
commodity 3 
F1d. 40 


iy the standard L.P. formulation, the constraint matrix 
of the whole problem is in block diagonal form, with the blocks 


corresponding to the constraint matrices of the subnetworks 
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in Fig. 48 and coupled by rows expressing the capacity 
constraints. 

It 1s well known for the single-commodity transship- 
ment problem with equality of total demand to total supply 
that any one of the node flow conservation equations (cor- 
meomonding to rows in the L.P. constraint matrix) is linearly 
@ependent on the other conservation equations. This is equiv; 
alent to saying that by conservation of flow over the whole 
network, given supplies and demands at all the other nodes, 
Emeenet flow at 5 node i must equal . Witheadwe. regard t© 
eaom>. For a sink node i, this means that given demand and 
impr y im the rest of the network, the net amount of commodity 
G@emee1ted at the node is predetermined to be b.. 

Hie above fine Of “reasenine 1S true ton seach eomeume 
icone cworks in Fig. 48. 

Suppose, therefore, in each of those subnetworks, we 
Eneose node B* to be the node with the linearly dependent 
eaeatron. If we re-combine Ai, iB back into a single neace. 
with demand Ey bs) , flow conservation within each sub- 
network will compel an amount bs to be deposited at B, which 
imemeprecisely what the original problem requires. fhe result 
ema Fig. 49, which is a single-commodity transshipment net- 
work with joint capacities over appropriate arcs incident to 


fee ihe capacities satisfy the conditions in paragraph 9 of 


mime previous Section. 
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2k 
B ny (-bp) 
k=] 





Key for arcs 
with joint capacities: 


mie ee 
2 nT 
—H- Yap 
Fig. 49 
Fig. 49 may then be reduced to the standard capacitated 
tuansshiipment network in Fig. 50, using Transformation I]. The 
Griginal flows are directly identifiable and no change to the 
objective function is needed. If necessary, other Transfor- 
mations may be applied to uncapacitate it or convert it to a 


transportation network. 
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Fig. 50 


Note that the above construction is valid, independent 
Ga: 


a. Whether the subnetworks (and hence the original problem) 
feed transshipment or tranSportation structure. 


b. The presence of multiple arcs between any node and node B. 


€. Whether the subnetworks are identical or not (except of 
Course for nodes Bk), This allows permissible arcs to be 
different for each commodity (provided, of course, the per- 
missible arcs for different commodities which are included 
in a joint capacity between a node pair are oriented in the 
same direction). 


d. The type of nodes (source, sink, or transshipment) in- 
volved. 


e. The number of commodities. 


imme Existence of lower bounds on arcs. If present, these 
gold be removed at the transformation stage depicted in 
Fig. 49, by using Transformation VIII. 


wee ine presence of uncapacitated arcs incident with B. 
minese would be dealt with by uncapacitated arcs directly 
between node B and appropriate nodes in each subnetwork in 
tie. 49. 
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oe, Complex Joint Capacities 
iemceseripeten jJOint'’ in the multicommodity context 


can be taken to apply in two dimensions, because:- 


a. For a given arc, the capacity need not include the sum 
of all commodities flowing along that arc. 


b. For a given commodity, the capacity may be imposed over 
Mmeet OL arcs ({1.€., jointly, in the single-commodity sense). 


If we imagine the subnetworks in Fig. 49 to be stacked 
Meena pile, then the concept of a two-dimensional joint 
Capacity leads to an interesting visual aid which can be used 


in small problems to check whether complex joint capacities 


at the busy node violate conditions for unimodularity. 
commodities 





arcs 
(RE) (2-By sees) 


Fig. 51 PIG aay 


Suppose Fig. 51 represents the plan view of the stack 
of subnetworks, and imagine all the arcs corresponding to 
Gliese incident #vith node B to be cut at plane P as shown. A 
cross-section in the plane would appear as in Fig. 52, where 
the circles represent cut arcs converging, say, into the page 


Gomme incident with node B, and the crosses represent arcs 


49 





Oriented in the opposite direction. Non-permissible arcs 
for any commodity would not be represented. 

Any complex joint capacity over arcs and commodities 
may be represented simply by sets of appropriate markings 
in Fig. 52. Such a diagram may be called a "constraint dia- 
gram" for the busy node, and may be examined for violation 
mee Conditions set out in Section II.C. Existence of an 
equivalent network model is not possible, if at any time: 


a. A set contains markings indicating arcs oriented in 
Soposite directions, or 


b. Two sets overlap (1.e., have an intersection which is 
a proper subset of both sets). 


Pie. Soe rliidsetrates both violations described above. 
Whee Fig. 54 illustrates constraints which make a network 
formulation possible. The geometrical outlines in the figures 


enclose arc sets which are included in the respective joint 





Sepacteies. 

a, N capacity 1 wv, capacity 1 
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arcs 
Fig. 53 Fig. 54 
hiershould pbeamoted a8aimvthnat ares Fepresented wan 


these diagrams are not node-disjoint, since they are all 
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ieiadent with B. Capacities therefore will not violate the 
Benadition III of Section II.C. 

For large problems, of course, a systematic computer- 
executed examination of capacities is more appropriate than 
such diagrams. 

Joint capacities which satisfy conditions for uni- 
modularity in the constraint diagram will also satisfy them 
when the problem is depicted as a single-commodity problem, 
oom Fig. 49. Transformations may then be used to simplify 
Someuex jOiInt Capacities aS required. 

feextensions of the Result 

If a multicommodity network has more than one busy 
modem(t.c¢., has some capacitated arcs which are node-disjoint) , 
Mempnomot interest to determine conditions which allow trans- 
formations to an equivalent unimodular single-commodity network. 

It is useful to review briefly why the addition of a 
second busy node would in general make an equivalent single- 
commodity transformation impossible (if this were possible, 
the general multicommodity problem would be unimodular). 
momsmaexs Fic, 49 jn the constructive proof of Result I. The 


basic fact is that the single-commodity formulation makes no 


k 
B 


total demand E(-b§) at node B of a single hypothesized com- 


distinction between the demands b It recognizes only a 
modity. What makes the busy node problem work is that flow 
conservation within the kt) subnetwork forces the Dreeuse 
amount bk to be deposited at B. This feature is in general 


B 
destroyed by the addition of another busy node, say, A. The kth 


D1 





subnetwork will then by flow conservation deposit a total 


of (ba a by) at nodes A and B, but not necessarily bs 


node A and bs at node B. An attempt to ensure that precisely 


by patil be deposited at node B by imposing a joint lower bound 


at 


Ot br en the algebraic«sum of flows along@the set of ares 
between node B and the kth network, for every k, will violate 
Condition I set out in Section II.C. Thais will be 
Seemeammediately from a constraint diagram for node B. 

In looking for extensions, therefore, it is necessary 
to look for conditions where: 


a. The capacities are such that they may be dealt with 
entirely within a commodity subnetwork. 


b. Busy nodes can be "collapsed" into an equivalent single 
busy node by eliminating arcs, usually at the cost of in- 
creasing the complexity of joint capacities and alterations 
in unit costs along remaining arcs. An appropriate tool 
for this purpose would be the Transformations (or their 
maverses) in the previous Section. A simple example is 
maeven below, while Result II derived later is a more complex 
@ase Of ‘'collapse."' 

meme velivery of the precise amount of each good can be 
ensured at each busy node, again usually at the cost of 
more complex capacitating arrangements. 

In attempting (b) or (c) above, it is again necessary 
Peomereek that capacities which may be added in the attempted 
Baichtormation dO not violate the three Conditions, either 
between themselves or with existing capacities in the original 
Pimowl em . 

An example of (a) above is in the following two 
closely-related special cases, where transformations to uni- 


modular equivalents exist for capacitated arcs at nodes other 


maome cin dddition to) the busy node. These cases arise when: 
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a. Each capacity at a node is specified to apply only to 
particular commodity. The constraint diagram for such 
a node could appear as in Fig. 55. 


seb) 


b. Joint capacities are applied across all commodities, 
mtepermissible arcs for each Capacity exist only for one 
commodity, as in Fig. 56. 
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Poem Capacity, such as the above at node(s) other than the 
busy node, are simply dealt with using Transformations, entirely 
Within the commodity subnetwork to which it applies. 

A simple example of "node collapse," is as follows. 


Consider the structure in Fig. 57 below with two busy nodes 


A and B, and r = 2 commodities, say. 
OWT Up3 AY 
CL 2 
ae ee VAB 
eo 
A TAB? TAB 
< Goh A 
Gok U 
4B 
Oy Yen w 
PiGeo! 


The structure may be transformed into the equivalent 


multicommodity structure in Fig. 58, with capacity Uap added 


ay 





memeres (1,;A) and (2Z,A) jointly, whose commodity unit costs 


k 
have CAR added on as shown. ay oe + at set have to be added 


to the objective function to bring in the cost of sending 


supplies at A to B, and flows i. are obtained from xk + 


1A 
me, Ak 
2 a 
U 
1A 
Cat od U3 AD 
6 U , . 
iA va AB an 
ait EO! 
+ 
2 + Ope U 
c AB 
(ayy, A 4) 


Fig. 58 

Fig. 58 is now a single busy node problem with a nested 
capacitating arrangement which can be transformed using Trans- 
formation V, after the problem has been broken up into the 
usual single-commodity subnetworks. 

In transformations of the above kind, it is necessary 
to see that conditions exist which ensure that the flows x\, 
in the eliminated arc AB are non-negative in the correct 


@azection. For the above example, if either: 


are reversed, it would yield, 


cpmeeearc (1,A) and/or arc, (2,A) 
K = k + ay (if both are reversed), 


me example, XB “Xny 7 X 


k 
A2 


Node,A became a source node, this would yield a = 
xk =x - pk A 
Lek A2 A 


In neither case would non-negative of a be assured after 


are AB is removed by node collapse. Note, however, that 


k 
Mmewersal of arc (B,4) would ensure non-negative Xap regardless 


: ; k b k 
_ = = ae 
Pimmome Orientation at A, since in that case Xap Xp3 Xp 4 
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If there are demands and/or supplies at both nodes A 
and B in the original problem, then their alhcebrale Sums ac 
node A in the transformed problem. 

An illustration of conditions which ensure delivery 
of the required amount of commodity at more than one busy 
node is as follows. Fig. 59 shows two busy nodes , Avand 5. 


part of a multicommodity problem with r = 2, say. Other parts 





Fig. 59 


Fig. 60 shows the problem broken up into commodity 
subnetworks, with A and B external to the subnetworks, and 


Transformation II applied to the Capacitated ares. 


Ist subnetwork > oi Ss Uaq) als 


Fig. 60 
; k k 
In Fig. 60, the kth subnetwork will Jd€cCept a total or a, + ap 


but not necessarily in the correct proportions from A and B. 
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To overcome the problem, joint lower bounds are placed on 


Ppposite arcs, as in Fig. 61, to ensure the proper amounts 


are supplied to nodes is and is 


cca 
| 






Fig. 61 


Transformation VIII can then be applied to remove the 
lower bounds. Note that this works only because each busy 
mea@emsupplies the rest of the network through only one node. 
I wthas were not the case, the lower bounds would be applied 
across node-disjoint arcs, which violates Condition II in 
feeeron [1.C. In this problem, either or both A and B could 


be sinks, and there may be more than two such nodes. 


C. RESULT II--THE TWO-SINK (OR SOURCE) TRANSPORTATION PROBLEM 
Evans, Jarvis, and Duke [4] have shown that an r-commodity 

(r > 2) transportation problem with 2 Sinks (or 2 sourees)) sand 

with all arcs individually capacitated, can be formulated as 

an uncapacitated single-commodity transportation problem. 

Figs. 62 and 63 show the transformation for a problem with 


3 sources, 2 sinks and 2 commodities. 
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Fig. 62. The Original problem, with 
2 commodities and 2 Sinks, 
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Fig. 63. The re-formulated, 
single-commodity problem. 
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In this subsection, the original multicommodity trans- 
portation problem is transformed into a capacitated single- 
commodity transshipment problem, which is considerably simpler 
than the single-commodity transportation formulation in Fig. 
63 above. It may be used in this form with suitable network 
eodes, or it may be further transformed into an uncapacitated 
maemsportation problem identical to that in Fig. 653. 

ert ine with the original problem (Fig. 62), construct 
Yr separate subnetworks, one for each commodity. Application 
of Transformation VII to each and every arc in the subnetworks 


mauewas Fig. 64. 


Ol Fe 
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Fig. 64 
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Note that: 


Peet sl inethe Onieanal problen ee them the Capacity 


X35 Sa; would be made redundant by the capacity Xo S Use, 

5 to be Dinding weer 

for r 2 2, both capacities are necessary, ee < - 
iG 

each subnetwork, and zx < Uss applied across all sub- 


Bence wWwevalready have U.. < a. for U- 
15 it of 
within 


2 nena | 
networks. 


RY 


ak k 
i wathy bounded capacity Key 4 Us, Se 
k 
Bele x2 ,; and hence is sequal cies co ce < U 
k=] 1 i Te “I =e 
which is required in the original problem. 


= 


14’ 

Ss. the particular requirement for two vsinks (oueEwo 
sources} in the original problem makes transformation to 
hoe. 64 possible. 

Fig. 64 is now equivalent to a single busy node problem 
with the sink node designated as "busy," with the dependent 
mass flow conservation constraint equation in each subnetwork. 
An analogous procedure to the "busy node" proof results in 
Fig. 65, a single-commodity upper-and-lower bounded trans- 
smipment problem. Costs along arcs incident with node Bare 


ZerO, oubSCripts refer to the original problem. 
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To get rid of the lower bounds, use Transformation VIII 
Oieene appropriate arcs to obtain Fig. 66, a single-commodity 
Capacitated transshipment problem with demand nodes at A, B, 
C, and D. Since we generally are not interested in the mag- 
nitude of resulting flows along arcs incident with node B 
and since costs are zero along these arcs, Fig. 66 may be 
iieoerly used 1n suitable capacitated codes without further 
post-optimal adjustments to flows along these arcs, or to the 

2 


a6 
Megeective function (apart from adding yz Z ave as stated 
k=] i=l ~ ? 


fipeeeic. 64). 
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Fig. 66. Single-conmodity capacitated transshipment problem. 


Finally, to obtain the uncapacitated transportation form, 
apply Transformation IV to each capacitated arc, to obtain 


ries O/ . 







@ Vis” (% ay-Uy 4) 


k 
LG 4) 
1.2 
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Pld en.6/ 


On re-arranging Fig. 67 (which has a transportation struc- 


fure) into visually bipartite form, we shall obtain precisely 
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2 
5) 


by equality of total demand and supply. 


k _ - 1 Kees 
Pave. 65 with es U.4) (b, +b by zUsy at node B, 


. 

It should be noted again that: 

ee ihe Cransshipment problem in Fig. 66 has significantly 
fewer nodes than either the original multicommodity problen, 
Oieeme Original end result in Fig. 63. Since solution times 
@re primarily related to the total number of nodes in a 
Meoohem, it can be expected that Fig. 66 with a suitable 
Capacitated code should result in considerable savings in 
ine. 

fine result 15 really a special case of the previously 
freeved extensions to the original "busy node" result described 


ijectime previous subsection D. 
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IV. UNCAPACITATING A TRANSSHIPMENT 
NETWORK--COMPUTER IMPLEMENTATION 
A. GENERAL CONSIDERATIONS 

ite cificiency of nétwork-based codes lies in Carrying out 
Miemateps Of the SIMPLEX algorithm by manipulation of spanning 
[eeec>s representing successive basic solutions. At the heart 
Oretnmese codes are efficient means of using list structures 
Mmemoreore and change a tree structure in computer memory, and 
to store and manipulate data required by the SIMPLEX algorithm 
nein the tree structure. 

Mie Transformations in Section I] change the original ner. 
perkeby the addition or deletion of nodes and arcs, and by 
emamees in costs and capacities of arcs. In a network-based 
Gage, corresponding changes are necessary to the list struc- 
tures which define the network and contain the data needed for 
the SIMPLEX algorithm. 

Contemporary efficient network codes are very fast, and 
so the area of interest lies in large problems where computer 
time and storage requirements are significant. The most 
efficient codes store three numbers for each arc, one of which 
emits capacity. Since most real problems have many more arcs 
than nodes, there is a considerable saving in computer storage 
her Uncapacitated problems. 

Interest is therefore drawn to Transformation IV, which 
mieapacitates an arc at the cost of adding two additional nodes 


and arcs. This obviously makes it unattractive for heavily 


63 





Capacitated problems. Its advantage is likely to lie in 
problems with only a small percentage of capacitated arcs, 
mere not having to reserve storage for all the are capacities 
more than offsets the slight increase in size of the overall 
network due to application of the Transformation. There may 
be a penalty, hopefully small, in speed due to the work needed 
to actually implement the Transformation in the code. However, 
1f a version of the code is written to handle uncapacitated 
Eweplems, then this time penalty is likely to be more than 
offset (depending on how lightly capacitated the problem is) 
fapene reduction in time due to simplified solution logic and 


data manipulation. 


peer eClS OF THE IMPLEMENTATION 

Im this Section, network problems with appropriate char- 
acteristics are chosen, and the Transformation implemented by 
Wepeaply modifying an existing computer code. Particular 
moecers Of the investigation are: 

1. Mechanical implementation of the Transformation, l.e., 
Omieesthe Original problem is input to the code. Detection of 
conditions in the problem which require applications of the 
Transformation, and the consequent applications, are done auto- 
iemeecally in the code. This applies also to the extraction of 
optimal flows and objective function of the original problem 
moon the solution to the uncapacitated problem. 

2. Where in the code to apply the Transformation. There 
are a number of places in the code where the changes of list 


structures due to the Transformation can be introduced. 
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Seeewie effect on the solution time caused by application 
Sietie Transformation and the variations in subparagraph 2 
above, compared with the case when the Transformation is not 
eppiied. It as to be hoped that speed improves, or at least 
not significantly decreased, 

meenough the results of the investigation hold only for 
the particular combination of Transformation, problem type, 
am code used, there is sufficient similarity with other 
Transformations or other codes to make the approach used in 


fice investigation to have quite general application. 


MeoeocCRIPTION OF THE CODE 

GNET [10] is a very fast primal network-based code in 
FORTRAN, designed for large-scale problems by Bradley, Brown, 
and Graves, which solves the general capacitated transshipment 
Beeptlem (and consequently, the other more specific network 
models) using the bounded SIMPLEX algorithm. 

l\ecy solution technique iS built on the structure eo. pace 
Polmeiions, which are known to form trees of basic arcs (ignor- 
ing orientation) spanning the nodes. The primal pivotal 
transformations of bases are performed by re-arranging, adding 
to, and deleting from, these trees structurally, rather than 
by numerical operations on coefficient matrices. 

A detailed description of its workings and its design 
philosophy is found in Bradley, Brown, and Graves [ |e 
brief description of its lists or "arrays" 1s given below, 
but this Section is written under the assumption that the 


reader has [2 ] available for reference. Note that the ole 
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variable names do not necessarily agree with notation used in 


mrevious Sections. 
Mapenre Data Structure 
This group of arrays defines the network structure of 

the problem. Each arc and its attached information is defined 
beeen entry into every one of these arrays, and they constitute 
the required input data to GNET. Supplies and demands are 
mionmeusing “dummy"’’ arcs connected to a fictitious super supply 
and super demand node respectively, but these are discarded 


after use in the preprocessor subroutine. 


M An integer defining the number of nodes in 
the problem. 


TC ) Aneare-length array Of staal! Rodeo (nc as 
nodes which arcs are directed away from). 


Ft) A node-length array of head nodes. 
Gl ) An arc-length array of unit costs of flows. 
Crt ) An arc-length array of are capaci t lesen 


a €Apacitated problem, arcs May sbeustomeanone 
of the basis with flow equal to Capacity, 
such arcs being marked by a negative sign 
(-CP( )) for identi treation: 


The reason why H( ) is only of node length is because 
all arcs with the same head are stored in contiguous space. 
The arc list may be viewed as sublists of arcs, each sublist 
belonging implicitly to a particular head node, followed by 
another sublist belonging to the next larger head node number. 
The jth entry in H( ), gives the address in the are list of 


tie first arc with head j. 
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UC ) 


P( ) 


SI eLex Data Structure 


(iemeconstraints of the dual problem to the primal 
transshipment problem are of the form: 


CG) ¢- -UC)) U1) = 0 for each basic -arem( ia 


The dual variables U( ), also called SIMPLEX multi- 
miverms, are stored in a node-Leneth array.) Sebeerune 
Mimare With tavorable SIMPLEX cost criteria to eneen 
PiemcUnrent basis involves, in GNEI, scaleulaeane fic 
iaietactor (abso called “2. - G2 “Singh. P ieemminiede cy) 
(imac ie) = Uy) + U7) forJeach arc in a candidate 
ist. and choosing that arc with the mose megariyceure 


It ethe tree representing a basis is Viewedoas “hanging 
dewnetrom a “root'’ node (the name for suchrarstrmecure 
iomarponescence' in graph theory). then a umeane 
"backpath'' exists from every node back to the root. 
laws, every node (except the root) has a singlespmuece 
decessor node, defined as the first node encountered 
OM this backpath. The minimum information needeaere 
mepresent a tree 1S a node-length array Py) )oreciee 
ment predecessors. The orientation Of an arc @iieeic 
tree 1s indicated by the sien of its predecessor. 
=P(i}) indacating an arc (1,P(1)), and +P( i), andwes. 
ome (P(i),i). .An important use of the predeces an 
list is in updating the flows in anew basis. If 

flow 1S changed in the entering are (1,)), thepem 
aeemangve of flow only in the arcs of the basis wenn 
are in the unique path from 1 to } im the eurpene 
tree. This path 1S idéntifred by jtterating vee 
predecessor relationship along the backpaths of i 

and j until a single node ‘is common to both. 


iMremcurrent tlow X( ject teaden basic arc —(i) eee 
aacociated with (i.e... sndexed by) the deepensnode 
of the pair (see explanation of array D( ) below). 
For example, flow along basic arc (i,j) is stored 
ime X( i) at jy 2S the predecessor of 1. 


For purposes of flow updating, an efficient method 
Shadentifying tie patho eroemanode 17 to sy s1neacic 
ourment tree. wiere (iejjeisetnc 2reuvomem cc irc 
Gurrent basis, is to store for each node the number 
of modes On the current backpath to the root; vealiied 
the depth D( ) of the node in the current tree. 
Depth information indicates which of i,j 1s deeper 
and should be iterated using P( ) along its back- 
path. When both backpaths are at equal demi, seme 
nodes are compared for identity. A match indicates 
the join node has been reached, otherwise both back- 
paths are iterated for another comparison. D( ) 1s 
also of use in identifying all the nodes belonging 
fowdny particular Suptice. | 
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e( ) 


et ) 


supp OSse —albasis arc (1,3) leaves the Current pasis. 
Witememj 15 the deeper nodevor the palms) selon. 

as part of the pivot sequence, only those U( ) 
associated with nodes in the entire subtree rooted 
Aegeave updated rom "the cCUurmeme pasi ss LEdeneinode 

ime tme entire tree. is given a Current Successor node 
number IT( ), which enables nodes in any subtree to 

be identified by exhaustively traversing through all 
its nodes starting from the root. The successor 
Retattonship is recursive in fiat, for any moce aes 
MeGuUroive SuUCCeSSOrs are Exhaustively traversed Demome 
any other nodes of the entire tree, and the rule 
adopted here is that leftward recursive successors 

are traversed first. A traversal (also called "pre- 
Order") through any subtree in Fig. 68 clarifies this. 


CPX{ ) is a node-length array containing the capacity 
minus the current flow of each current basis arc. 
SUppese arc (1,)}) 1s to enter the basis. ee laemeeee 
ouikgomme arc needs to be identified as therarc sco. 
meavondime €O the Minimum of : 


eee es) 


be ~X(€_) alone the bDackpath fromm togene 
Tea node 


c. CPX( ) along the backpath from j to the 
FOI node. 


ihis is the ‘ratio ™tesw wot the SEMPLEX -alecomernne 


Figure 68 below illustrates the main SIMPLEX data 


fatayevariables in relation to a basis arborescence. 








@2) x (12) 
CPX(12) 
@) No Fig. 68. Continued next page. 
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Node i: 
P(1) 
D(i) 
Tete ) 
(9,2) 
GPX (3): 
X (3) 
X(12) 
CPX(12): 


Backpath from node 9: 9, 4, 8, 


Backpath from node 2: 2, 7, 


12 if 6 5 


Incoming arc 


6 i 8 
f= Selah 
4 2 Z 
9 i 4 


Capacity minus flow; for are (5411) 


Flew allome arc (53511) 


Plow aiearc (€2.12) 


Capacity minus flow, for amem@ 12) 


ai 


ee 


Unique path from 9 

to 2-18 9; 43a 
iy 2 With Le elie seroma 
node. 


eee’ ,11) is outgoing, U( ) for only sub-tree 7, 2, 12, 10 


are changed. 


and sets up an initial arborescence which consists of 
ficial arcs between each node and the root. 


of these arcs is set equal to the demand or supply of 


GNET imtroduces an additional mode ‘caliecdmene 


root, 
arti- 
The flow in each 


the 


associated node (arcs are directed to the root from sources 


and from the root to sinks). 


wane appears typically as in Fig. 


OT; 


The initimal arborescence stsuc- 








The rest of the SIMPLEX data structure is set up, with sinks 
given an arbitrarily large negative value (“Bie M") fore team 
Bie). Sources have zero U( ). 

pee eceprocessor LNET 

GNET includes a preprocessor subroutine LNET which 
Nas as its main purposes: 

a. To check for certain input errors, and produce 
Statistics on the problem, such as total demand and supply, 
number of nodes and arcs, ranges of costs and capacities, 
percent of arcs capacitated (if any), type of problem (capaci- 
tated/uncapacitated, transshipment/transportation), etc. 

Bee lO Set Up the arc data structure im headenoede 
sequence as described previously. 

Gay 1O attach a X()) valuevto cach node caua meromine 
Supply Or demand, which are used by GNET as initial basis 
flows, as described above. xX( ) at sinks are marked with 
negative signs, which are subsequently removed in GNET after 


their use for identification purposes. 


DD: THE TRANSFORMATION CHOSEN FOR IMPLEMENTATIOW 
Transformation IV was chosen for the reasons in subsection 
A above. Transformation IV uncapacitates an arc, as shown, 
iepome rig. 70 to Fig. 71, by adding two artificial nodes and 
meres. Arc (K,j) is chosen to be regarded as containing the 
Oepeinal ('real') flow, although (1,K+1) may also be used for 


this purpose. 
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F low=X F low=X 
Cost=C C Cost=0) 


GD K+} Demand= -CP 


F low=CP-X 
Cost=0 


CK) Supply=cp 
F low=X 
Cost=0 
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aS THE CHOSEN PROBLEMS 

For reasons given in subsection A above, this Transforma- 
tion is senn to be most efficient in lightly capacitated 
problems. 

There are some well-known large test problems produced 
by a program called NETGEN and reported in Ref. 5, which are 
widely used to test codes. Two of the largest capacitated 
problems are in fact lightly capacitated and are suitable for 
this investigation. Their main characteristics are tabled 


below. 


yA 





NETGEN 39 NETGEN 40 
Total demand (= total supply) 4,000,000 2,000,000 
Number of arcs 15,000 25,000 
Number of nodes 5,000 3,000 
BemGentage of arcs capacitated O97 es 0.66% 
Number of pure sources/sinks 1se7 700 100/300 
Number of transshipment sources/sinks 79/300 49/100 
Maximum incoming/outgoing arcs to Aes 207 26 
a node 
Cost range MIG J EtG, 
Capacity range 5, 0010 ZOOS 

Oe Ors 5 Zaire 


Rabie I. TEST PROBLEM CHARAGIERFS@ ses 


These large problems require a lot of computer memory space, 
SO program modifications were tried out on a small problem 
with 7 nodes before implementation on the large problems. 

If the storage space for arc capacities need not be used, 
eimseresults in respective savings of 15K and 23K single pre- 
cision numbers, corresponding in the IBM 360 to 60K and 92K 


bytes respectively, which is substantial. 


PF, BASIC DATA CHANGES 
ie Transformation introduces changes to the are data 
Seructure, from Table II to Table III, which correspond to 


ieee 70 and Fig. 71 respectively. 
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AEE, Last 


Location Index 2 T(%) C(%) CP(%) 
HS) 
ue ; phe a T(L)=i c(L) cP (L) 
H(j+1)-1 


NAR = Total number 
Onmares 


Table II. ARC DATA STRUCTURE--UNTRANSFORMED PROBLEM 








/NG eae BRS o 
Mocation Index ¢ 1aE>) C(2) CP(2) 
1 
6 am wo ; : 
; Er T(L)=K CCL MAXCP 
(K,j) oi : 
NAR 
: {Arcs 
added 
by 
A , . ; : Transtor- 
eapeen fw sttons 
Arc 28 ; ’ MAXCP 
(i, ea Oe ale ; | | 


Table III. ARC DATA STRUCTURE--TRANSFORMED PROBLEM 
is 





Two additional entries into the arc list are made for each 
Capacitated arc (i,j) as shown, while the original arc (av) 
becomes (K,j). All arcs are uncapacitated by setting CP( ) = 
MAXCP, an arbitrarily large number > total supply. H(K) = 
H(K+1), since K is a pure source node, i.e., its arc sublist 


meneyacuous. 


G. INSERTION OF DATA CHANGES INTO THE PROGRAM 

There were three places in the program where the Transfor- 
mation was experimentally implemented. 

1. Within LNET, where every capacitated arc is trans- 
formed. GNET consequently accepts and treats the added arcs 
and nodes as if they were no different from those in the un- 
transformed problem. 

2. AS in 1 above, but with modification to GNET so that 
the initial arborescence is set up with each Transformed arc 
Memceaing a structure as in Fig. 71, instead of having all 
nodes (including the added nodes) at the same depth as normal 
meal (Fig. 69). This is done for all capacitated ares 
before pricing-out for the first pivot, and may be termed 
Swme-hanging.” 

3. “On-the-fly" in GNET, where the Transformation is 
spplied to a capacitated arc only when it is about to enter 
ampasis. 

Since the problem is uncapacitated in all the above cases, 
there cannot exist non-basic arcs with flow at capacity, 1.e., 
CP( ) is never marked with a negative sign. Two program 


statements which tested for this condition were removed from 
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the pricing-out procedure. This was the only modification 
to solution logic which took into account the uncapacitated 
nature of the Transformed problem. 
A description of each of the three cases follows. 
me implementation*®Within LNET 
This was the simplest application. After the arc 
data for all the arcs in the capacitated problem had been 
structured as normally in LNET, each arc was tested to see 
mel was Capacitated and if so, the data changes in Table [11 
were inserted. Each source node added by a Transformation 
had positive X({ }) value attached, equal to the capacity of 
the Transformed arc. Sinks had negative X( ) values. 
GNET does not differentiate between original and added 


nodes, and sets up the initial arborescence normally as in 


fee / 2. 
_ Root Mt2J+1 
J = number of 
capacitated 
arcs. 
Last added 
Added nodes node M+t2J 
ilgwe2e 
2. "“Pre-hanging" 





In this case, changes to data structure occur in LNET 
exactly as in 1 above. In addition, within the main program, 
the initial arborescence is set up so that the two nodes k and 
k+] which were added in LNET by Transforming capacitated arc 


(i,j) hang down from the original tail node 1. This is done 


AS 





for every capacitated arc, and amounts to re-structuring nodes 
Mei to M+2J in Fig. 72 above to hang below the appropriate 
Original nodes. A typical initial arborescence will look as 
in Fig. 73 below, where node 2 originally had two Capacitated 
outgoing arcs, for example, and nodes 3 and M are not the tail 


nodes of any capacitated arcs. 


Oo” & ©: -O..- 7 
D> © © © 
D> O° © 


Fig. 73 

Ii the initial arborescence was as in Fige 7/2 seme 
Memere a Gapacitated are ({1,;})) canventer a basis, two prior 
Pers tO bring nodes k, k+l into the tree must be performed. 
@eme-hanging'' as in Fig. 73 in effect saves the work needed 
to perform these two prior pivots, and it was thought that 
some saving in time due to reduction in total number of pivots 
to optimality could thereby be saved. 

GNET is allowed to set up the normal uniform-depth 
initial arborescence up to the Mth node. All arcs in loca- 
tions increasing from 2 = H(M+l1) are known to have been added 
ie contiguous pairs in LNET, one pair for each are uncapaci- 


@atred, as in Table III. Therefore, the following procedure 
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momadopted for each &, where & begins at 2 = H(M+1) and in- 
Mero s Dy Steps Of 2, (1.e., 2%, 2842, 244, ...)4 The first 
foae, kK Of an added node pair is found from k = T(2). The 
second added node is immediately obtained as k+l. The original 
tail node iis found from i = T (+1). All three nodes needed 
to pre-hang an added node pair are thus identified, and the 
Summeture is defined by setting the following structure 


Vewemaples (relate these to Fig. 73): 


P (k} Sorin. Ged 
P (k+1) =] 
eT k) = IT(i) If i1had more than 1 outgoing 


Capacitated arc; these two seate. 


ie (1) = k+] ments will cause subsequent node 
pairs to be pre-huns to) thewlews 
of preceding pairs: 

ae Oh ) = k 

D(k) = 3 

Dekel) = 2 

IT(M) = M+2J+1 (the root) 


inmeddaition, the following SIMPLEX data is set: 


CP X(k) = 0 

X(k+1) = -X(k+1), [Recall X( ) was set negative 
by LNET“tor-all sinks.4 

EPA(k+1) = 0 

Dek) = U(1) 

Wek) =U (i) 


The program then proceeds to the first pivot With tie leur weme 


basis corresponding to the pre-hung arborescence. 
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oS nem On- the= tly” 
This case is the logical extension of the two previous 
cases. LNET is modified only in order to mark Capacitated 
arcs by a negative T( ). The main program starts off nor- 


mally with only the original nodes and arcs (including capaci- 


tated arcs). The arc (i,j) which is chosen to enter a basis 
iemecsted for a negative T( ). If it is a capacitated arc, 
then: 


a. The two additional nodes and arcs are created and 
eorresponding changes made in the arc data structure as in 
ase | above. 

b. The two additional nodes, k, k+l, are pre-hung 
meee the arborescence as in case 2, from node i. The are to 
fimeer Che basis is now (k,j). 

€. Appropriate SIMPLEX data is attached touemewodged 
nodes. 

The above procedure is essentially a simultaneous 
combination of cases 1 and 2, performed between the pricing- 
out and pivot sequences in the main program. The motivation 
behind hang-on-the-fly is to minimize the number of added 
nodes in bases by creating them only when a capacitated arc 
is about to enter. This should reduce the overall amount of 
work done in arborescence traversing, and updating the arbor- 
eseence and SIMPLEX data in the pivot sequence. The size of 
the network and hence the basis arborescence grows as the 
number of pivots increases, since nodes added by transforming 


entering capacitated arcs remain as part of the network. 
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This is in contrast to cases 1 and 2, where the network is 


expanded to full (constant) size at the outset. 


fee RESULTS AND DISCUSSION 
ime ecrtect of the Transformation on the speed of solution 


mietne various cases examined is shown below. 


NETGEN 39 NETGEN 40 
Time Number of Time Number of 

(seconds) Pivots (seconds) Pivots 
Unmodified GNET Tua age U555 G70 6409 
ferpeacitated problem) 
Transformation within iors 9615 7058 (SII 
eee 
Pre-hanging PAs) 9689 Gi 2 Oe 
Hang-on-the-fly ILLES Re: S19) 3) 4 64.7 O7 55 


Table IV. COMPUTER IMPLEMENTATION RESULTS 


GNET is written in standard ANSI FORTRAN. All computer 
runs were on the Naval Postgraduate School's IBM 360/67. The 
solution times given below are actual execution times, not 
Cock times. 

Meris clear that the results are quite unitorm across aim 
cases in speed of solution or number of pivots in the above 
meses, for each problem. It is possible to discern a siiugnt 
improvement by hanging-on-the-fly compared to the other cases. 

Thus, speed performance is maintained even with the Trans- 
formation implemented. This is significant in two respects: 

1. If a version of GNET is developed to specifically 


solve uncapacitated problems, solution speed sno) a be “at 


fies: 








least maintained while there would be a considerable saving 
miesteorage presently required to store arc capacities. This 
would enable still larger or more heavily capacitated problems 
to be accommodated. 

2. Similarly, an uncapacitated GNET version would mean 
Miemetron Of the solution logic and data handling designed 
specifically for the capacitated problem. This should give 
a speed improvement compared to the same problem solved in 
@im@eecurrent GNET. 

Wyerali, therefore, an uncapacitated GNET should solve 
Barger Or more heavily capacitated problems for little speed 
loss, or the same (lightly capacitated) problem at a greater 
speed. 

Of course, if the code available can only solve the 
Miteapacitated transportation problem, then use of this Trans- 
formation becomes a matter of necessity rather than efficiency. 

Por progressively more heavily cCapacitateds roblens. ene 
saving in storage requirements due to deletion of CP({ ) in an 
uncapacitated version of GNET would be increasingly offset by 
the storage needed by the expanded size of the problems. Each 
application of the Transformation adds two more nodes and arcs. 
In the extreme case of a 100% capacitated problem, the number 
m@eearcs will expand by a factor of three, and this is elearly 
well beyond the point where application of the Transformation 


would be considered worthwhile. 
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